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Abstract 

We show that the two-component system of hyperbolic conserva- 
tion laws dtp + d x (pu) = = dtu + d x p appears naturally in the 
formally computed hydrodynamic limit of some randomly growing in- 
terface models, and we study some properties of this system. 
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1 Introduction 

The macroscopic behaviour of physical systems can often be described in 
terms of non-linear partial differential equations. In many cases, it had 
been shown that functionals of microscopic models from statistical physics 
converge in the hydrodynamic limit towards certain solutions of these partial 
differential equations. 

Studying the partial differential equation (or the system of partial dif- 
ferential equations) can turn out to be a very hard challenge in itself: Ap- 
pearance of singularities in finite time, shocks etc. The so-called hyperbolic 
conservation laws have in particular received a lot of interest. Even in one 
space dimension, these PDEs proved to be extremely interesting and chal- 
lenging both mathematically and phenomenologically. These are partial 
differential equations of the type 

dtu + d x J{u) = 

where u = u(t, x) takes its value in W 1 and J is a non-linear function from 
W 1 into R n . 

The best known and most investigated examples are the following. (See 
e.g. {?], [yj, 15] for a comprehensive introduction and survey of the subject.) 



(1) Burgers' equation (with no viscosity): n = 1 and 

d t u + d x (u 2 /2) = 0. 



(2) The isentropic gas dynamics equation in one space dimension: n = 2, 
the components are the density field p(t, x) and momentum field m(t, x) 

| d t p + d x m = 

[ d t m + d x (m 2 / p + p{pj) =0 

where is the pressure, depending on density only. 

(3) The so-called p-system, which is an alternative formulation of the dy- 
namics of one-dimensional gas. The two components are the velocity 
field u(t,x) and the specific volume (= inverse density) field v(t,x): 

{ dtv — d x u = 

{ (2) 
[ d t u + d x p(v) = 0. 

Here p(v) denotes the pressure, as a function of specific volume. 

(4) The shallow water equation is another two component system: h(t, x) 
denotes the height of the (shallow) layer of water, u(t, x) is the velocity 
field: 

f d t h + d x {hu) = 

\ d t u + d x {u 2 /2 + ti) =0. 

Since Riemann, a considerable amount of knowledge and technology 
(more recently, for instance, entropy solutions, compensated compactness 
method) has been derived that give a better understanding of the physically 
relevant solutions to these equations. 

In the present paper, we will be considering a particular two-component 
(i.e. n = 2) system of hyperbolic conservation laws that arises in the context 
of surface growth (or more precisely growing interfaces, since the surface is 
one-dimensional). In other words, at each time t > 0, one sees a landscape 
x i— > h(t, x) where i£i The function h is increasing in time. The rough 
phenomenological description of the phenomena we are interested in corre- 
sponds to the case where the surface is growing in the normal direction to 
its boundary, but there exists a 'tension' that tends to keep the surface to- 
gether, in the sense that it will fill in holes quickly. In the physics literature, 
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a famous equation has been proposed by Kardar, Parisi and Zhang (the KPZ 
equation) as a model for such situations, cf. ||. It is (in the mathematical 
jargon) an ill-posed non-linear partial differential equation with a stochastic 
term: 

d t h = Ah- (d x hf + W 

where W = W(t, x) denotes a space-time white noise. We do not want 
to give a review of the huge physics literature on this equation, but we 
briefly stress two aspects. (See Q] for a state-of-the-art survey of the physics 
literature on the subject and an exhaustive list of references up to 1995.) 
First, there exists to our knowledge no completely satisfactory (see however 
||) derivation of this equation from a microscopic model. Second, it is 
predicted that 'the' solution to this equation has a special scaling behaviour 
at late times. More precisely, it is believed that when a, t, x are very large, 
h,( a \t,x) = a~ l ^h{at, a 2//3 x) is also a solution to the KPZ equation. The 
exponents 1/3 and 2/3 should be related to various conjectures and recent 
rigorous results concerning the fluctuations of highest eigenvalues of random 
matrices, of first passage percolation paths, of longest increasing sequences 
etc etc. 

One way to define one-dimensional interfaces hit, x) in terms of particle 
systems goes as follows: Start with a (finite or infinite) system of particles 
that evolve randomly in the potential h(t, x) (or in some potential defined 
in terms of h) and that all contribute to increase the potential in the sense 
that h(t, x) corresponds to the joint local time (i.e. cummulated occupation 
time density) of the particles at time t and site x. In other words, h(t, x) 
increases locally at x if there is a particle at x and time t. Note that 
this leads naturally to a two-component system in the (formally computed) 
hydrodynamical limit: the first component is the density of particles, and 
the second component is the gradient of the profile of the potential. 



In [17], we constructed a continuous stochastic process, corresponding 
on a heuristic level to the case where there is exactly (and only) one particle 
(its location at time t is denoted by Xt) which is driven by 

dX t = -d x h(t,X t )dt 

and h(t, x) is the local time of X at x and time t, so that 

d t h(t,x) = 5(X t -x). 
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For details concerning the construction and primary properties of this pro- 



cess and a rigorous version of these equations, see |17|]. Let us just empha- 
size a couple of features: The process (Xf,t > 0) is a random process, even 
though the previous 'differential equations' look very deterministic. One rea- 
son is that (in the stationary regime), the function x i— > h(t, x) is not regular; 
in fact, it is a Brownian motion in the space variable (for fixed t). Second, 
Xt is not a usual stochastic process (it is not solution of a stochastic differ- 
ential equation for instance), it has the 2/3 scaling: (a~ 2 / 3 X at ,t > 0) has 
the same law as (Xt,t > 0). In particular, (a _1//3 /i(a 2//3 x, at), t >0,i£l) 
has the same law as (h(x,t),t > 0, x £ M) so that h has the same scaling 
property as the asymptotic scaling conjectured for the KPZ equation. 

The process (Xt,t > 0) can be viewed as the scaling limit of a discrete 
negatively reinforced (i.e. self-repellent) random walk (S n ,n > 0) on Z 
called the 'true self-avoiding walk' in the physics literature. This is a nearest- 
neighbour walk on Z that decides at each step to jump to the left or to the 
right according to a probability depending on how many times it has visited 
the neighbouring sites (or edges) before. Suppose for instance that after n 
steps S n = x and that the discrete walk (Si)i< n has jumped already I (resp. 
r) times on the edge immediately to the left (resp. to the right) of x. Then, 
S n +i = x + 1 with probability 



e & 



P(5 n +i = x + l\l,r,S n = x) = - 



e~ 



-0r 



where j3 > is some fixed constant. In other words, the walk will prefer 
to go along the edge it has visited less often in the past. Note also that 
the probability in fact depends only on the difference I — r (which depends 
on all the past trajectory). The distribution of the rescaled position of the 
random walker, 5 n /n 2//3 , converges to (a multiple of) the one-dimensional 
marginal distribution of the continuous process X t described above, fl6( |. 

It seems natural to consider the case where this one particle is replaced 
by many particles performing the same kind of self-repelling walk on Z, with 
a joint cumulated local time of all particles. Or, in the continuous space-time 
setting: a continuously distributed cloud of particles (that all contribute to 
the same local time), which is the subject of the present paper. As we 
shall see, this leads in the (formally computed) hydrodynamic limit to the 
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following system of hyperbolic conservation laws: 

d t p + d x {pu) = 

( 4 ) 

d t u + d x p = 

where p corresponds to the density of particles at x and time t, and u(x, t) = 
—d x h corresponds to the negative gradient of the interface. It seems, that 
although this system looks very natural, it has not been considered in the 
literature. We should emphasize that in spite of some formal similarities with 
the p-system (Q) and the shallow water equation (||), the system (H) shows 
very different behaviour and describes a quite different phenomenon. We 
hope that its study may lead to improved understanding of some aspects of 
'growing interfaces' in general. In particular, this equation could shed some 
light on some of the conjectured properties of the KPZ equation. The goal 
of the present paper is not to present a complete treatment of this system 
of partial differential equation, but rather to initiate it as an alternative 
approach to 1-d domain growth and deposition phenomena. 

2 The PDE: phenomenological derivation 

We define a deposition model in the following terms. The actual state of 
the system is described by two functions: 

p : R + x R -> R + and h : R + x R -> R. 

p(t, x) is the density of the population performing the deposition, while 
h(t, x) is the deposition height at time t and space coordinate x. The rules 
governing the time evolution of the system are the following 

(1) The total population is conserved, so that the continuity equation 

d t p + d x (pu) = 

is valid, where u(t, x) is the velocity field, to be specified by the dynam- 
ical rules. 

(2) The deposition rate is proportional to the density of the population, i.e. 

d t h = cip, (5) 
where c\ is a positive constant. 
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(3) The population is driven by a velocity field proportional to the negative 
gradient of height 

u = -c 2 d x h, (6) 

where c 2 is another positive constant. This rule corresponds to the 
self-repellence mechanism described in the introductory section. 

From (H) and @ we readily get 

d t u + cic 2 d x p = 

Without loss of generality, we can choose c\c 2 = 1 and get the two compo- 
nent system of hyperbolic conservation laws 

dtp + d x (jm) = 
d t u + d x p = 

This system of PDEs with initial conditions 

p(0,x) = p (0) {x), u(0,x) =u {0 \x) (8) 

is the main object of the present paper. 

As a first remark we mention here the scale invariance of (j^). Let v € K 
be fixed. Given the functions (t, x) i— > p(t, x) and (t, a?) i— > n(t, x) and a 
positive fixed number q, define the rescaled functions 

p (a) {t,x) :=a 2 ^p(at,a v x), 
u (a) (t,x) := a l - u u{at,a u x). 

One can easily check that if (p,u) is solution of ([?]), then (p^ a \ mW) is also 
a solution, for any a > 0. The choice v = 1 yields the hyperbolic scale 
invariance valid for any hyperbolic conservation law. More interesting is 
for our purposes the choice v = 2/3. This is the physically relevant scale 
invariance, since the density changes covariantly under this scaling, i.e., the 
total mass J p( a 'dx is unchanged. 

With this choice of v the following scale invariance of the deposition 
height follows: 

h^(t,x) :=a- 1/3 h{at,a 2/3 x). 

Recall that this is exactly the conjectured asymptotic scale invariance of the 
one- dimensional KPZ equation. 
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3 Bricklayers 



We define a system of interacting particles living on Z, with two conserved 
quantities, whose hydrodynamic modes are governed by a two-component 
system of hyperbolic conservation laws which, after taking another limit (low 
density/late time), transforms into our system (Q). The computations of the 
present section are somewhat formal. Working out all technical details (e.g. 
proving uniqueness of the equilibrium Gibbs measures or technical details 
of Yau's hydrodynamic limit) needs more effort. The present section serves 
as microscopic motivation of the PDE proposed above. 

3.1 The particle system 

The Great Wall of China is being built by a brigade of bricklayers. The wall 
consists of columns of unit-size bricks, piled above the edges of the lattice Z. 
The height of the column piled above the edge (J, (i.e., number of bricks 
in this column) is hj. In the dynamics of the system the discrete negative 
gradients Zj := hj—i — hj £ Z will be relevant. The bricklayers occupy the 
sites of the lattice. At each site j E Z there might be an unlimited number 
rij € N of bricklayers. Bricklayers jump to neighbouring sites and at each 
jump j — > j ± 1 a brick is added to the respective column of bricks. 

In more technical terms: particles (= bricklayers) perform continuous 
time nearest neighbour walk on the lattice Z and hj measures the cumulated 
(discrete) local time on the lattice edge + 1). 

About the dynamics: the jump rates are chosen so that the following 
conditions hold: 

(1) the bricklayers' jumps are driven by the local shape of the wall so that 
they try to reduce the differences Zj (i.e. to keep the height of the wall 
even) , 

(2) conditionally on the actual shape of the wall the bricklayers jump inde- 
pendently. 

This is done as follows. The instantaneous rate of jump from site j to 
site j ± 1 (for each bricklayer sitting at site j) is equal to r(±Zj), where 
r : Z — > (0, oo) is a fixed monotone increasing function which defines the 
model. In order to be able to compute explicitely the stationary measures 
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3 bricklayers at site j 

n =3 
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each of these 3 bricklayers 

jumps to the left with rate: r(4)/ V(~4) 
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j-2 j-1 j j+1 j+2 
Figure 1: The great wall being built 



(see subsection 3.2) we impose that r(l — z)r(z) is a positive constant (this is 
for instance the case if r(z) = exp(/3z)), and multiplying time by a constant 
term, we can in fact restrict ourselves to the case where 

r{z)r{-z + 1) = 1, for all z € Z. (9) 

Thus, the following changes of configuration may occur: 

(n-j, Zj), (nj+i, Zj+i) -> (n-j - 1, Zj - 1), (n j+1 + 1, z j+1 + 1) 
with rate rijr(zj), and 

(n-j, Zj ), (nj-uzj-!) -> (n-j - 1, Zj + 1), (n^-i + 1, - 1) 

with rate rijr(—Zj). 

Clearly, Ylj n j an£ i Ylj z j are formally conserved quantities of the dy- 
namics. It is also clear that besides these globally conserved quantities the 
parity of rij + Zj is also conserved on each lattice site jeZ. 
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Now we give a more formal description of our interacting particle system. 
For s G {0, 1} let 

(N x Z) s := {(n, z)eNxZ:n + z = s mod 2}. 

Given the sequence s = (sj)jgz G {0, 1} Z we define the state space of our 
system as 

a, :=n( NxZ )./ 

Elements of fi s will be denoted by cj, i.e. to = (ujj)j^z with ujj = (rij,Zj) G 
(N x Z) g . The (formal) infinitisimal generator of the Markov process de- 
scribed verbally in the first paragraph of this section, is: 

L/H = J2n j r(z j )(f(e j+ u) - f(u)) + J2n j r(-z j )(f(O j ^) - f(u)), 



where the maps Qj + and 0j_ act on the subsets {w e Sl s : rij > 1} as 



respectively 



( (ni,Zi) 

(rii -l,Zi- 1) 

, (Ui + l,Zi + 1) 



(ni,Zi) 

(rii -l,Zi + 1) 
, (rii + l,Zi- 1) 



if 
if 
if 

if 
if 
if 



i = j 
i = j + 1 

« 7^ 3,3 ~ 1 
t = j 
i = j - 1 



3.2 Equilibrium Gibbs measures 

For k > denote 



i2(z) :=[I^) 



fe=i 



and 



0* := lim r(k) G (l,oo]. 
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Note that (||) implies that for all z E Z, 

R(-z) = R(z) = R(z - l)r(z) = + l)r(-z). (10) 

Fix the parameters s £ {0, 1}, A > 0, G (1/6*, 6*) and define the probabil- 
ity measure /i Sj A,6» on (NxZ) s as follows: 

, 1 A n 0* 

Z S (A, 0) n! /t(z) 

where 

A n 0* 



Z B (A,0):= £ 



n! 

is the normalizing factor (partition function). The measure fJ. s ,\,8 is a prod- 
uct measure on N x 7L restricted to the subset n + z = s mod 2. It is worth 
noting that 

Z s (\,6) = Z s (A,0 -1 ) and fJ, s ,\,o(n, z) = li s ,\,e-^(n, -*)■ (H) 

For a fixed sequence s G {0, 1} Z and fixed parameters A > 0, G 
(1/0*, 0*) we define on O s the probability measure 



By direct computations, one can check using (10) that for any function 
/ that depends only on the value of finitely many uj^'s, for any fixed Sj, Sj+i 
and 



n j r(z j )n S3tXfi (uj j )^ Sj+ltXfi (uj j+1 )f{Q j+ (u:)) 
^2 n j+1 r(z j+1 )fj, s . t x,e(^j)ti Sj+1 ,x,e(^j+i)f(^) 



and a similar dentity holds for the jumps to the left. It follows that given 
the local parities nj + Zj = Sj mod 2, the probability measures /x s \ : q are 
stationary for the dynamics. These are the equilibrium Gibbs measures of 
our system. For a similar computation in the context of a simpler one- 
component domain growth model see also Q . 

Invariance under spatial translations is unfortunately lost in this very 
general setup. In order to impose it, we restrict ourselves to one of the 
following two choices: either s = or s = 1. 
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3.3 The hydrodynamic equations 

For the rest of this section we fix either s = or s = 1 and we do not denote 
any more the dependence on s. 

As we have mentioned already the globally conserved quantities of our 
system are ^ ■ rij and Y^j Zj. In the equilibrium regime fix,e the- averages of 
these quantities are 

dfogz(\e) Q diogz(\,e) 

P ■= {nj)x,0 = A — , u := {zj)x,e ~ 




(12) 



dX ' Jl ' 89 

These are the particle density (per site) and the average slope of the height 
of the wall, in equilibrium. It is easy to see that the map K+ x (1/6*, 9*) 3 
(A, 9) i — > (p, u) G E + x R is globally invertible. Indeed, 

/ dp/dX dp/ 89 \ _( Var(n) Cov(n,z) 
I du/dX du/d6 J \ Cov(n,z) Var(z) 

So the gradient matrix on the left hand side of fll2| ) is everywhere invertible 
and this implies global invertibility of the map (X,8) i— > (p,u). With slight 
abuse of notation we denote the components of the inverse function A = 
X(p,u) and 9 = 9(p,u). From (|ll]) it follows that 

X(p, —u) = X(p,u) and 9(p, —u) = l/9(p,u). (13) 

In order to guess the system of hydrodynamic equations we have to see 
first how the infinitisimal generator acts on the conserved quantities. An 
easy computation shows: 

Lrij = (rij-ir(zj-i) - njr(-Zj)) - (rijr(zj) - n j+1 r(-z j+1 )) 
Lzj = (rij-ir(zj-i) + n j r(-z j )) - (rijr(zj) + n j+1 r(-z j+1 )) 

On the right hand side of these equations we see discrete gradients of fluxes. 
This fact helps us guessing the hydrodynamic equations. Applying the stan- 
dard formal manipulations to our gradient system (see e.g. ||, H) and using 
the straightforward identities 

(n i r(±z j ))A,e = A0 ±1 

in the hydrodynamic limit taken with hyperbolic (Eulerian) scaling of space 
and time, we arrive at the system of PDEs 

r d t p+d x (x(p, u )(9(p,u)-9(p, u r 1 )) = o, 

\ d tU + d x (x(p,u)(9(p, u ) + 9(p, u )- 1 )) =0. 1 
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Under growth conditions on the rate function r(z), as z — > oo, Yau's 
'relative entropy method' (see e.g. |lg], |§, ||]) in principle can be applied 
to our system of interacting particles, resulting in the validity of the above 
system of PDEs in the hydrodynamic limit, as long as the solutions are 
smooth. 

From the system (llj|) we can derive the system (|7j) by taking a second 
limit: We replace p(t, x) by a 2 / 3 p(at, o?l 3 x) and u(t, x) by a 1//3 ii(ai, a 2 / 3 x) 
We note that for small values of the variables p and u, 

\(p,u) = p + o(p), 0(p,u) = 1 + cu + o(u). 

where 



d 2 z 



d6 2 



G (0,oo) 



A=O,0=1 / 

Letting now a — > 0, we arrive at (]?]). We should emphasize here that this 
scaling limit does not depend much on the details of microscopic system. 
Also, from any conservation law of the form 

J d t p + d x J(p,u) = 
\ d t u + d x K{p,u) =0 
we would get (Q) under the same limiting procedure, provided that 

J(p, u) = pu + o(pu), K(p, u) = p + o(p), as p, u — > 0. 

This indicates that (0) is valid for a wider class of microscopic systems. 

4 Analysis of the PDE 

We are now going to see how the methods developed in the PDE literature 
(see 0, [l4|, 15]) can be applied to our system. In order to put things into 



perspective, we briefly recall general results and see how they can be applied 
in the context of our system (Q). 

4.1 Two-component systems of hyperbolic conservation laws 

For a generic two-component system we shall use the notation v = v(t, x) = 
(vi(t, x), V2(t, x)) T . (The superscript T will denote transposition of vec- 
tors/matrices.) The generic two-component system is 

dtv + d x J(v) = 0, (15) 
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where v i— > J(v) = (J\(v), J2{v)) T is a smooth vector field over RxR, J is 
the flux of the flow of the conserved quantity v. The initial conditions are 
specified by 

v(0,x) = v (0) (x). (16) 

For a (possibly vector- or matrix valued) function / = f(v) we denote 
the gradient with respect to the v- variables V/ = (df/dv±, df/dv2)- For 
classical smooth solutions v(t,x), ( Jl5| ) is equivalent to 

d t v + V J • d x v = (17) 

(we use • to indicate products of matrices). 

As a technical device one usually also considers the so-called viscous 
equations 

d t v + V J • d x v = ed 2 x v. (18) 
Existence and unicity of smooth solution v^ £ \t,x) of for any bounded 



and smooth initial conditions (16) is guaranteed by the smoothening effect of 
the artificial viscosity term on the right hand side. One hopes that physically 
acceptable (stable) solutions of the original system ( |l5|) can be obtained as a 
strong limit of the viscous solution v^ £ \t, x), as e — > 0. The existence of this 
strong limit is a very difficult problem and is a main object of investigation 
in the context of hyperbolic conservation laws. 

In our case (|7|) the two components are v\ = p, v% = u, and the corre- 
sponding fluxes are Ji(p,u) = pu, J2(p,u) = p. The inviscid system is (0). 
The (artificially) viscous system is 

f d t p + d x (pu) = ed 2 x p ^ 
[ dtu + d x p = ed x u. 

The viscous solutions (which do exist and are unique) will be denoted by 
(p< e )(t,s), ««(*,*)). 

4.2 Hyperbolicity 

One has to check that the matrix V J has two distinct real eigenvalues p < A. 
The domain where this holds will be denoted 

£>h yP := {v G M x R : p(v) < \{v)}. 
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The corresponding left (row) and right (column) eigenvectors will be 
denoted by I and r, respectively, m and s. That is: 



I ■ VJ = A/, V J • r = Ar, (20) 

m ■ VJ = pm, VJ • s = fj,s. (21) 



For our system we find: 

VJ 



u p 
1 



and 



A = +i( v W4p + u ), Z = (A,p), r = (A,l) T , (22) 

= --(yV + 4p-«), m=(n,p), s = (p,l) T . (23) 

Note that / • s = m • r = 0, as it should be. 
We conclude that for our system, 

D hyp = {(p, u) G R x R : u 2 + 4p > 0}. 

Note that in the physically relevant domain with non-negative densities 

V ph :={{p,u) eMxK:p>0}, 

there is one single point where strict hyperbolicity is lost, namely (p, u) = 
(0, 0). On the other hand, we found that the system is still hyperbolic in the 
physically meaningless domain X\yp \ Pph = {(p,u) € ^hyp : P < 0} ^ 0. 
At the moment nothing seems to prevent solutions to flow into this domain. 
Later we shall see that Lax's maximum principle (valid for stable entropy 
solutions) takes care of this problem. 

4.3 Riemann invariants, characteristics 

In the generic two-component case, we are looking for scalar functions 
2\yp 3 v ' ^ w(v) S R and space-time trajectories R + 3 t t— > £(t) G R 
such that for smooth solutions of ( |T5| ) (or, equivalently, of (|17|)) w is con- 
served along the trajectory £(t), i.e. 

^«>(u(t,£(t)) =0. 
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Using ([l?]) we find: 

d£ _ (Vro ■ V J) • d x v 
dt Vw ■ d T v 



(24) 



In order to solve (24), Vw must be a left eigenvector of the matrix VJ. 
It follows that this relation admits two solutions: one for each eigenvalue of 
VJ. We denote the two solutions by w (corresponding to the eigenvalue A), 
respectively, by z (corresponding to the eigenvalue p). The gradients Vw, 
respectively Vz, are parallel to the row vectors I, respectively m, defined in 
(|20|), respectively (^)- In other words, 

Vw-s = 0, ^ = A, 

Vz ■ r = 0, — = a. 

dt r 

These equations, of course, do not determine uniquely the functions w(v) and 
z(v). Given two smooth, monotone maps /, g : R — > R, the transformation 
w := f(w), i := g{z) leaves the above equations invariant. The functions 
w and z are called the Riemann invariants, or characteristic coordinates of 
the problem. 

In our case the most convenient choice of the Riemann invariants w and 
z is the following: let 



V w := {(p, u) G V hyp : ^u 2 + 4p - u > 0}, 
V z := {(p, u) € V hyp : ^u 2 + 4p + it > 0}, 

and define w : P„, — > R, z : P 2 — > R by the formulas: 



w(p, u) = — y \J u 2 + 4p — u ^\/n 2 + 4/> + 2u^j 



= - V v 7 ^ 2 + 4p + u \^/u 2 + Ap- 2u) . 



Note that T> w nT> z = T> p h, so that both Riemann invariants are defined 
in the physically relevant subdomain. 

It is straightforward to check that both Riemann invariants w and z 
defined above are convex functions of the variables (p,u). This fact will 
have crucial importance in later analysis. 
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Figure 2: Level lines of the Riemann invariants: z = est and w = est 
4.4 Genuine nonlinearity 

In plain words, genuine nonlinearity of a two-component system of hyper- 
bolic conservation laws means that on the level curves w(v) = const., re- 
spectively z(v) = const., the characteristic speed li, respectively A, varies 
strictly monotonically. Formally: 



dw 



oz 



Performing straightforward computations this turns out to be equivalent to 

VA • r ^ ^ Vll ■ s. 

That is: the characteristic speeds A and \x vary strictly monotonically in the 
direction of their corresponding right eigenvectors. 

In our case, given the formulas (p2[) and (p3[) we easily get 



m 2A „ 2u 

A — /i (J, — A 
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Recall from (^), ( p3[ ) that on D p h we have \x < < A, with strict inequalities 
for p > 0. We conclude that our system is genuinely nonlinear in the interior 
of the physically relevant domain 2? p h- On the half lines p = 0, u < 0, 
respectively, p = 0, u > (on the boundary of I'ph) genuine nonlinearity of 
the first, respectively, of the second, characteristic speed is lost. 

4.5 Weak solutions, shocks, Rankine-Hugoniot conditions 

As it is well-known, a nonlinear system of hyperbolic conservation laws (|i~5| ) 
can develop singularities (e.g. discontinuities), irrespectively of the smooth- 
ness of the initial conditions. A generalized or weak solution of (|15|), ( jl6| ) 
in a space-time domain is a bounded, measurable function (t,x) h- » v(t,x) 
satisfying 

/OO /"OO 
/ (a t 0(i,x) -v(t,x) + <9 x 0(t,x)-J(v(t,x))}cftdx + (25) 
-oo io 

/oo 
0(0, x) • t>W(x)dx = 
-oo 

for any row vector valued test function = (0i,02) with compact support 
in the respective space-time domain. This last equation is obtained by a 
formal integration by parts. It is easily seen that a strong (smooth) solution 
is also a weak solution. 

Assuming a (locally) piecewise C 1 solution with a spatially isolated jump 
discontinuity at some space-time position (t, x) £ R+ x R, one derives the 
Rankine-Hugoniot conditions which relate the left- and right limits of the 
function x \— > v(t,x) at the discontinuity and the propagation speed of the 
discontinuity: 

J 1 (v(t,x+))- Ji(u(t,x-)) J 2 (v(t,x + ))-J 2 (v(t,x-)) 

? 7 = O" = — — , (2b) 

Vl(t,X + )—Vl(t,X ) V2(t,X + )—V2(t,X ) 

where a is the propagation speed of the discontinuity, i.e. the slope in space- 
time of the line of discontinuity. (|2^) is derived from (^5|) by an elementary 
local argument, using the divergence theorem (in space-time). Given the 
two independent relations in (^), any three of the five values vi(t,x~), 
i>i(i,x + ), V2(t,x~), i>2(i,x + ), a determine the other two. This imposes a 
serious restriction on the possible jump discontinuities of weak solutions. 
Note that the conditions are left-right symmetric. 
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We turn now to our system (|7]). We denote by (p left , u left ), respectively 
(/9 nght , u nght ), the values of the component functions at the two sides of the 
presumed discontinuity. The Rankine-Hugoniot conditions are: 

„, right „right _ „, left Jeft ^right _ Jeft 

" r " H _ _ P_ P (27\ 

pright pLcit ^right ^left ' V ' 

Given the value at one side of the discontinuity, the value at the other side 
as function of propagation speed is expressed as follows: 

left 

p right = a 2 -au Mt , u ^ = a -t— (28) 

a 

Note that p nght , respectively, it right , is expressed as function of a and u Mt , 
respectively, as function of a and p left , only. (In principle, both should be 
expressed as functions of a, p loft and it left .) This is a special feature of our 
system. 

The propagation speed, as function of the values of the components on 
both sides of the discontinuity, is expressed as: 

a± = ±^{^/(u right ) 2 + 4p lcft ± u right } = ^{^/(u lcft ) 2 + V ight T u Mt }. 

Lax's condition of stability for Rankine-Hugoniot discontinuities, |lCt |, 
specified for two-component systems reads as follows: Assume that the weak 
solution ( |25| ) of the two-component system ( |l5| ) is piecewise smooth, with 
a spatially isolated discontinuity with values f; left , respectively, u right on the 
two sides, propagating according to the Rankine-Hugoniot conditions (|26|). 
The discontinuity is a stable back shock, respectively, front shock, according 
whether 

^(v right ) < a < mm{fi{v Mt ), A(v right )}, (29) 

or 

max{A(t/ ight ), p(v Mt )} < a < \(v Mt ). (30) 

Rankine-Hugoniot discontinuities which do not obey either one of the con- 
ditions ( p9[ ) or (|30|), are unstable, physically not realisable. 

Tedious (but, in principle straightforward) computations show, that in 
the case of our system (|j) the discontinuities propagating according to ( |27| ) , 
or equivalently (|28|) are stable back shocks if a < and stable front shocks 
if a > 0. 
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4.6 Entropies 

Given the two-component system of conservation laws (p"5|), we look for 
additional conserved quantities, i.e., for pairs of functions I'hyp 3 v *—> 
(S(v),F(v)) elxR which satisfy 

d t S(v) + 8 x F(v) = (31) 



for smooth solutions of the original problem fllq ) (or, equivalently: for 
smooth solutions of (O)). Indeed, ([H]) means, that S(v(t,x)) is globally 
conserved quantity, with flux F(v(t,x)). The pair of functions (S,F) is 
called entropy /flux pair. Using the form (|i~7|), valid for smooth solutions of 



(15), one finds the system of PDEs defining an entropy/flux pair: 

VF = VS • V J, (32) 

or, in extended form: 

dF_ = yds_dJi_ k = l2 

dv k dvi dv k ' 

This is a two-component linear hyperbolic system of PDEs for the two un- 
known functions S and F — just well determined. There are various alter- 
native equivalent ways of writing it. E.g., eliminating the function F we get 
a second order hyperbolic PDE (a wave equation with variable coefficients) 
for S: 

dJ 1 d 2 S fdJ 2 dJA d 2 S dJ 2 d 2 S 

+ ~ T: I „ „ + 7i o" = 0. 



dv2 dv\ \dv2 dvi J dv\dv2 dv\ dv\ 
Or, changing variables to the characteristic coordinates (w,z): 

dw dz ' dz dw 
Or, eliminating F between these two equations: 

d 2 S 1 /dfidS dXdS 



dwdz A — /x \dw dz dz dw 

These last two forms explicitly show the wave-character of the entropy equa- 
tions (|32"|), Of particular importance are those entropy /flux pairs for which 
the function v i— > S(v) is convex. Such pairs will be simply called (with 
slight abuse of terminology) convex entropy/flux pairs. 
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In the case of our system (^) the entropy equations, written in terms of 
the physical variables p and u, are: 



dF_ _ dS_ dS_ d£_ dS_ 
dp dp <9w ' du ^ dp 



Or, eliminating F: 



a 2 s u d 2 s d 2 s _ Q 3 

dp 2 dpdu du 2 

The existence of a strictly convex entropy/flux pair, globally defined on 
Dpii = {(p, u ) '■ P > 0, u € R} and with S bounded from below is very 
important, since the applicability of Lax's Maximum Principle cited in the 
next subsection relies on it. Here it is: 

u 2 

S(p,u) = plogp+ —, F(p,u) = up(\ogp+ 1). (34) 



Lax's 'entropy wave construction' (cf. applies also to our system 

(^). Since these computations are rather involved, we do not reproduce them 
here. Let us just point out, that this robust method ensures the existence of 
a sufficiently rich family of convex entropy/flux pairs in any fixed subdomain 
compactly contained in P p h- 

There are also other (more ad hoc) methods of constructing entropy/flux 
pairs. Following, e.g., the ideas of |l3| we may try to find so called similarity 
solutions of the entropy equation (p3|) of the form: 



S(p,u) = p a <f>{pPu) (35) 

Elementary computations show that (3 = —1/2 is the only choice consistent 
with (p3[). Inserting (p5[), with (3 = —1/2 into ( |33"[ ) we find the following 
ordinary differential equation for the function <fi : R — > M: 

3(y 2 - 4/3)^"(y) + (5 - 8a)y<f/(y) + 4a(a - l)4>(y) = 0. (36) 



Any solution of (pq), with any a £ R fixed provides an entropy of our 
system, via (|35|). So, we are able to construct a sufficiently rich family of 
entropy /flux pairs to our system (^). 
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4.7 Entropy solutions 

A weak solution (t,x) t— > v(t,x) of the generic system ( |T5| ) is called entropy 
solution if for any convex entropy /flux pair (S, F) we have 

d t S + d x F < (37) 

in the sense of distributions, i.e., for any positive test function (t, x) t— > </>(£, x) 

/OO /'OO 
/ {d t (j){t, x)S{v(t, x)) + ^0(t, x)F(v(t, x))} dtdx > 0. 
-oo JO 

Entropy solutions are the only physically admissible, stable ones among 
the weak solutions. Strong limits of all convergent approximation schemes 
(such as vanishing viscosity or various convergent finite difference schemes) 
result in entropy solutions. It is also expected that convergent hydrody- 
namic limits of interacting particle systems result in entropy solution of 
the corresponding hyperbolic conservation laws. For piecewise smooth weak 
solutions, Lax's stability condition for the shocks mentioned in a previous 
paragraph is equivalent with the entropy conditions (|37|). 

Of particular interest is the following Maximum Principle, due to P. Lax, 
see e.g. Q. 

Maximum Principle for Entropy Solutions. Assume that the follow- 
ing two conditions hold 

(i) The Riemann invariants v i— > w(v) and v i— > z(v) of the system of 
hyperbolic conservation laws fiH\ ) are (globally) convex functions of v. 

(ii) There exists a globally defined convex entropy /flux pair, with entropy 
function bounded from below. 

Then, starting with bounded initial data, sup_ 00<x<00 |u^°^(x)| < oo, along 
entropy solutions (t,x) i— > v(t,x) the maximum values of the Riemann in- 
variants, sup_ 00<x<00 w(v(t,x)) and sup _ 00<x<00 z(v(t, x)) do not increase 
with t. 

Remark: The same statement applies for solutions v^ £ \t, x) of the viscous 
system (|l8| ) — this follows from the classical maximum principle. If 
converges strongly as e — > 0, then the limiting v is in fact an entropy solution 
of the inviscid system ( |T5| ) and forcibly it obeys Lax's Maximum Principle. It 
is not clear whether all entropy solutions arise as limits of viscous solutions, 
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with vanishing viscosity. A general proof of the Maximum Principle for 



entropy solutions can be found in [11]. 



Applying this theorem to our system we find that if we start with 
bounded initial data x i— > (p^ (x) , u^°\x)) € 2? p h (that is: with non-negative 
initial density) then entropy solutions will stay in the physical domain, i.e., 
for any t > x i— > (p(t,x),u(t,x)) G D p h- (See Fig. 2 for graphical rep- 
resentation of the level curves w(p,u) = const, and z{p,u) = const.) This 
is a very important consequence of the Maximum Principle: as we already 
mentioned, a priori we could not see any reason banning a (physically rele- 
vant) solution from flowing out into the physically meaningless domain with 
p < 0. 

In the case of isentropic gas dynamics, (|l]), choosing convex versions of 
the Riemann invariants w and z, for any w max E M, z max £ the domains 

{(p,m) Gl+xl: w(p,m) < w max , z(p,m) < z max } 

are compact. So starting with bounded initial data global boundedness of 
(viscous and) entropy solutions is guaranteed by the Maximum Principle. 
This is unfortunately not the case for our system. The domains 

{(p,u) e T> ph : w(p,u) < w ma , x ,z(p,u) < z max } 

are not compact, see Fig. 2. So here is an open question: Is it the case, that 
if the initial data ^) are bounded then the solutions [p^ £ \t,x),u^ £ \t,x)) of 
the viscous equation fiH\ ) stay bounded for ever? Similarly: is it the case 
that entropy solutions of ^) with bounded initial data stay bounded? We 
guess that the answer to these questions are affirmative, but we could not 
prove this yet. 

4.8 Vanishing viscosity, existence of entropy solutions 

The existence of entropy solutions for a two-component syetem of hyperbolic 
conservationlaws (|l5[) is a notoriously difficult question. The most powerful 
approach seems to be the program initiated by R. DiPerna in Q , completed 
for the case of isentropic gas dynamics (pi) in H , then refined and extended 



in Lions et al. [13] and in several other papers. 



In @, DiPerna proves the following result: 
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DiPerna's Theorem. Consider the two- component system of hyperbolic 



conservation laws (15) and the corresponding viscous system (jig). Assume 
that 

(i) The Riemann invariants v ^ w(v) and v h- > z(v) are convex. (More 
precisely: there are convex choices of the Riemann invariants. See 
subsection 4-3.) 

(ii) The system is genuinely nonlinear. (See subsection 4-4-) 

Let C be a domain compactly contained in £>h yp and assume that the sequence 
of solutions v( £ \t, x), t £ [0,T], x G K, e — > 0, of the viscous systems 
ftl(\ ) takes values from C. Then there is a subsequence v^ £ >(t,x) which 
converges strongly in L} oc ([0,T] x R). The limit v(t, x) is entropy solution 
of the system (15). 
Some Remarks: 

(1) The proof relies on the construction of Lax's 'entropy waves', hinted at in 
subsection 4.6 and essentially on the so-called compensated compactness 
method developed by Murat and Tartar. We do not have a chance to 
reproduce here any technical part of the proof. 

(2) It is assumed that the viscous solutions stay in the domain C. However, 
even in this form the theorem is technically very- very difficult. Extra 
difficulties arise by relaxing this condition and imposing conditions only 
on the initial data: in the isentropic gas dynamics and in our case too, 
the solution data will typically flow to the boundary of the domain of 
hyperbolicity and genuine nonlinearity, p = 0, where this theorem is not 
any more valid. 

(3) For extensions, physically more satisfactory formulations and enormous 
further technical difficulties see e.g. j|, [13], etc. 

This theorem can be applied in a straightforward way for domains C, 
compactly contained in P p h- We can add to this that if initially 

maxw < or maxz 

X X 

then, due to the Maximum Principle, the viscous solutions p^ e \u^ are kept 
away from the 'dangerous' vacuum line p = 0, see Fig. 2. So, in this case 
one has to care only about the boundedness of the solutions. 
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